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ABSTRACT 


Proportional  Navigation  and  Coinmand  To  Line  Of  Sight 
missile  guidance  are  explored.  A  system  flow  graph  is 
developed  for  each  guidance  technique.  The  block  transfer 
functions  are  developed  and  a  state  space  representation  of 
the  systems  is  defined.  The  missile  systems  are  then  tested 
using  one  two-dimensional  engagement  and  two  three-dimensional 
engagement  scenarios.  The  final  three-dimensional  scenario 
introduces  measurement  noise  in  order  to  evaluate  the  effect 
of  noise  on  the  guidance  algorithms.  The  engagement  results 
are  then  compared  to  analyze  the  miss  distance  of  each  type  of 
missile  guidance. 
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I. 


INTRODUCTION 


A  guided  missile  can  be  controlled  using  two  different 
methods.  The  first  is  when  the  missile  contains  its  own 
guidance  system.  This  type  of  missile  is  beneficial  in  that 
once  fired  it  will  track  its  target.  The  second  type  of 
guidance  has  a  ground  fire  control  system  to  command  the 
missile.  This  type  of  missile,  called  command  guided,  does 
not  contain  a  target  seeker.  Two  radars,  or  one  radar  capable 
of  tracking  two  targets,  are  required  at  the  fire  control 
station;  one  will  track  the  missile  and  the  other  the  target. 
The  fire  control  system  will  calculate  the  required  missile 
acceleration  commands  and  relay  them  to  the  missile  by  either 
a  radio  link  or  fiber  optic  cable. 

The  type  of  guidance  system  implemented  is  largely 
dependent  on  the  missile's  mission.  A  long  range  missile  will 
need  a  self  contained  guidance  system.  A  point  defense 
missile  will  use  a  self  contained  seeker  or  command  guidance. 

The  guidance  system  supplies  the  input  to  the  missile 
control  system.  We  will  use  a  roll  stabilized  "skid- to-turn" 
missile.  The  roll  stabilization  will  permit  a  simpler 
analysis  because  there  is  no  coupling  between  pitch  and  yaw. 
Figure  1  shows  a  block  diagram  for  a  missile  control  system. 
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Figure  1 .  General  Missile  Guidance  System 


Missile  interception  simulations  using  command  to  line 
of  sight  and  proportional  navigation  guidance  systems  are 
developed.  Chapter  II  explains  the  guidance  laws.  Chapter 
III  develops  the  simulation  algorithms.  Chapter  IV  tests  the 
algorithms  with  known  two-dimensional  results  and  a  three- 
dimensional  problem  with  and  without  measurement  noise. 
Chapter  V  discusses  the  simulation,  conclusions,  and 
recommendations . 
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II.  MISSILE  GUIDANCE  LAWS 


A.  GENERAL 

The  missile  guidance  system  provides  the  autopilot  with 
the  necessary  information  to  produce  the  required  acceleration 
commands.  The  missile/target  intercept  geometry  has  several 
important  parameters.  Figure  2  depicts  a  typical 
missile/target  intercept  scenario. 


y 


Figure  2 .  Missile  And  Target  Intercept  Scenario 
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Several  important  parameters  can  be  developed  by 
analyzing  Figure  2 . 

Rn,  :  Tracker  to  missile  range 
Rt  :  Tracker  to  target  range 

:  Tracker  to  missile  line - of  -  sight  angle 
Ot  :  Tracker  to  target  line - of  -  sight  angle 
Ym  :  Missile  velocity  vector  angle 
Yt  :  Target  velocity  vector  angle 
The  two  guidance  techniques  to  be  discussed  are  proportional 
navigation  and  command-to-line-of -sight . 

B .  PROPORTIONAL  NAVIGATION 

Proportional  navigation  missiles  are  guided  by  either 
reflected  or  emmitted  energy  from  the  target.  A  passive 
missile  will  be  guided  from  the  IR,  EO,  or  RF  emmitted  by  the 
target.  An  active  missile  will  send  an  RF  signal  out  to  track 
the  target.  In  each  case  the  energy  is  received  by  a  seeker 
which  tracks  the  target . 

Proportional  navigation  is  based  on  the  rate  of  change  of 
the  missile  to  target  line-of -sight  (LOS) .  The  missile 
commanded  acceleration  is  proportional  to  the  rate  of  change 
of  the  LOS.  The  ratio  of  the  missile  turning  rate  to  the  LOS 
rate  of  change  is  called  the  proportional  navigation  constant 
(N)  . 


Y 

N'  =  ^ 

6 

Y  (2.1) 

N  ‘  —  N' 

V 

m 

The  proportional  navigation  constant  must  be  greater  than  2  to 
ensure  system  stability.  A  proper  value  of  N  will  ensure  that 
the  missile  to  target  angle  will  remain  constant  thus 
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ensuring  missile  intercept.  Figure  3  illustates  this  point. 


Assuming  an  acceleration  is  applied  at  right  angles  to 
the  velocity  vector  of  the  missile  for  a  period  of  time  dt, 
the  missile’s  velocity  will  then  be  V„(t+dt).  The  velocity 
vector  will  have  changed  direction  by  dVn,- 
Assuming  a  small  angle  approximation  yields 

a  dt  V  dsf 

m  m  *  m 

(2.2) 

a  = 

Combining  (2.1)  and  (2.2)  leads  to 

a  =  vm  (2.3) 

m  in 

This  result  is  the  proportional  navigation  law  that  will  be 
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implemented  in  this  simulation. 


C .  COMMAND  GUIDANCE 

The  Command  To  Line  Of  Sight  (CLOS)  missile  is  given 
guidance  commands  that  keep  the  missile  in  the  LOS  between 
the  launch  point  and  the  target.  The  distance  between  the 
missile  and  the  LOS  is  defined  as  the  cross  range  error 
(CRE) .  The  fire  control  system  will  supply  the  proper 
commanded  acceleration  to  drive  the  CRE  to  zero. 

Since  two  separate  radars  are  required  for  this  type  of 
guidance  the  problem  geometry  is  slightly  different  than 
previously  described.  Figure  4  shows  the  CLOS  system 
geometry . 


z 


Figure  4 .  Command  To  Line  Of  Sight  Geometry 
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III.  SYSTEM  DEVELOPMENT 


A.  OVERVIEW 

The  system  block  diagram  is  shovm  in  Figure  1 .  The  block 
transfer  functions,  system  dynamics,  and  simulation  equations 
will  be  developed  for  the  simulation. 

B.  RADAR  DEVELOPMENT 

Target  flight  is  tracked  using  angles  in  the  pitch  and 
yaw  planes.  The  pitch  plane  is  defined  as  the  vertical  plane 
that  contains  the  target  and  the  radar .  The  yaw  plane  is 
defined  as  the  xy  plane. 

1.  Proportional  Navigation 

Proportional  navigation  system  geometry  is  shown  in 
Figure  5 . 
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From  Figure  5  the  following  angles  can  be  defined: 


^myaw 

Missile 

yaw  angle 

n 

^mpitch  • 

Missile 

pitch  angle 

n 

^  tyaw 

Target 

yaw  angle 

^tpitch 

Target 

pitch  angle 

The  system  requires  that  the  following  ranges  be  defined: 

R„  :  Radar  to  missile  range 
Rt  :  Radar  to  target  range 
R  :  Missile  to  target  range 
By  applying  elementary  trigonometry  to  the  cartesian 
system  geometry  defined  in  Figure  5,  the  following  equations 
can  be  derived 
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The  radar  system  will  produce  the  following  angles 

Oyaw  ’  Missile  to  target  yaw  plane  angle 

o^picch  •  Missile  to  target  pitch  plane  angle 

* 

The  angles  are  given  by  the  equations 


a  =  arctan 

yaw 


V 


pitch 


arctan 


(3.2) 


The  radar  will  send  these  angles  to  the  respective  yaw  and 
pitch  seeker  elements. 

2 .  Command  Guidance 

The  CLOS  radar  will  produce  a  cross  range  error  signal 
and  relay  this  signal  to  the  missile  autopilot.  The  cross 
range  error  is  the  distance  between  the  missile  and  the  radar 
to  target  LOS.  Figure  4  shows  the  CLOS  geometry. 

From  Figure  4  and  vector  calculus  the  cross  range  error 
of  the  missile  can  be  defined  as  follows 


|C^|  = 


1^  X 
«  /n _ t 

i^.i 


(3.3) 


This  calculation  yields  the  following  equation 


\CRE\ 


l-^tl 


(3.4) 


The  missile  autopilot  requires  that  the  cross  range  error 
be  broken  into  the  yaw  and  pitch  components.  Analyzing  Figure 
4  yields  the  following  equations 
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(3.5) 


CRE  =  Jx  ^+y  ^  sin  (o,  -  a  ) 

■  ,/CRE^-CRE^.^ 

The  sign  function  ensures  that  the  pitch  plane  cross  range 
eziox  can  be  positive  or  negative.  A  positive  cross  range 
error  indicates  that  the  missile  is  leading  the  LOS .  A 
negative  cross  range  error  indicates  the  missile  is  trailing 
the  LOS. 

C.  SEEKER  DEVELOPMENT 

1 .  Proportional  Navigation 

The  seeker  for  proportional  navigation  measures  the  rate 
of  change  of  the  missile  to  target  LOS  angle.  A  simple 
gimballed  seeker  will  use  the  angular  rate  of  change  of  the 
seeker  head  as  an  estimate  of  the  rate  of  change  of  the  LOS 
angle .  Figure  6  shows  the  seeker . 
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The  equation  of  motion  of  the  seeker  head  will  yield  the 
estimate  of  the  angular  rate  of  change  of  the  LOS .  The  seeker 
head  equations  for  the  pitch  and  yaw  planes  will  be  identical . 
We  will  develop  the  equations  for  only  the  yaw  plane.  The 
equation  of  motion  for  the  seeker  head  is  given  by 


T  = 


I 


SBBker 


(3.6) 


where 

T  =  Torque  applied  to  the  seeker  head 
I  =  Moment  of  inertia  of  the  seeker  head 
P  =  Seeker  bore  sight  angle 
Solving  (3.6)  yields 

$■  =  Z  =  -kjO-o)  -  (3.7) 

Taking  the  Laplace  transform  of  (3.7)  gives 

s^3(s)  =  (s) (3.8) 


Then  we  solve  for  the  seeker  transfer  function 


3(s) 

0{s) 


(3.9) 


where  Tsh  is  the  seeker  head  time  constant. 

A  typical  seeker  head  time  constant  is  1/8,  using 

this  value  produces  the  following  constants 
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=  64 

16 


(3.10) 


The  signal  flow  graph,  using  these  constants,  can  be  seen 
in  Figure  7 . 


64  1/8  l/s 


-1 


Figure  7.  Proportional  Navigation  Seeker  SFG 


The  following  state  space  representation  can  then  be 
implemented 


X 


SH 


0 

1 

0 

-64 

-16 

SH 

64 

X 


SH 


J 

1 

2 

L 

u 


SH 


(3.11) 


U  =  G 

SH 
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2 .  Command  Guldeuice 

The  CLOS  missile  control  system  does  not  contain  a  seeker 
head.  All  missile  control  functions  are  processed  and 
developed  by  the  fire  control  system  located  at  the  radar 
site . 

D .  GUIDANCE  DEVELOPMENT 

1 .  Proportional  Navigation 

The  missile  guidance  system  implements  the  proportional 
navigation  law  explained  in  Chapter  II.  The  major  difference 
is  that  an  estimate  of  the  angular  LOS  rate  is  used  vice  a 
measurement  of  the  actual  LOS  rate.  Therefore,  the  rate  of 
change  of  the  missile's  velocity  vector  is  given  by 

y  =  N  $  (3.12) 

*  jn  ^ 


This  leads  to  the  following  state  variable  representation 


■ 

- 

■ 

’  m  pj.  tcti 

N 

0 

ft 

pitch 

y 

•  ill  5'aw 

0 

N 

^yaw 

(3.13) 


2 .  Command  Guidance 

The  guidance  for  a  CLOS  missile  is  developed  from  the 
rate  of  change  of  the  missile’s  cross  range  error.  The 
missile  acceleration  is  equal  to  the  rate  of  change  of  the 
cross  range  error .  This  rate  of  change  is  then  used  as  a 
commanded  acceleration  in  the  autopilot . 

The  commanded  acceleration  is  developed  to  provide  good 
niissile  response.  To  ensure  good  response  the  missile 
acceleration  must  be  of  the  form 

+  {OL  +  |3)s  +  cx3  (3.14) 
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This  will  provide  the  damping  necessary  for  the  missile  to 
perform  correctly. 

Using  equation  (3.14)  the  following  commanded 
acceleration  is  developed 


(3.15) 


(3.16) 


a  =  40  and  p  =  196  produces  two  real  roots  at  s=-5.7l7l  and 
s=-34 .2829  . 

The  signal  flow  graph  for  the  guidance  system  is  shown  in 
Figure  8 . 


Figure  8.  Command  To  Line  Of  Sight  Guidance  SFG 


A  state  space  representation  of  the  guidance  system  is 
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(3.17) 


a _ _ 

^^pitcb 

196 

40 

0 

0 

a _ _ 

dncly,^^ 

0 

0 

196 

40 

CRE  ^ 

pitch 

CJ^  . 

pitch 


CRE 


yaw 


CRE 


yaw 


E .  AUTOPILOT  DEVELOPMENT 

1 .  Proportional  Navigation 

A  simple  autopilot  can  be  developed  by  applying  a  torque 
about  the  center  of  gravity  of  the  missile.  Analyzing  the 
equation  of  motion 


T  = 


m 


(3.18) 


and  noting  that  this  must  also  satisfy  equation  (3.14)  to 
achieve  stable  performance,  yields  the  following  relationship 

V„  ■  ^  -  -kV.  •  KWg  (3.19) 

■^OS 


Taking  the  Laplace  transform  of  (3.19)  yields 

_  KN 

$ (s)  s  *  K 


(3.20) 


and  defining  T;^p  as  the  autopilot  time  constant  produces 


K  = 


(3.21) 


The  signal  flow  graph  for  the  autopilot,  with  k=l,  is  shown  in 
Figure  9 . 


15 


Figure  9.  Proportional  Navigation  Autopilot  SFG 


The  state  space  representation  can  be  written  as  follows 


-1  0 

1  0 

Y 

0  -1 

0  1 

L  J 

u 


AP 


1 

■ 

N  0 

$ 

0  N 

ft 

u 


AP 


(3.22) 


The  missile  acceleration  commands  can  be  derived  by 
looking  at  the  missile's  velocity  vectors.  Figure  10  shows 
the  two-dimensional  missile  acceleration  geometry. 
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Figure  10.  Missile  Acceleration  Geometry 


It  can  be  shown  that  the  velocity  in  the  pitch  and  yaw 
planes  is  given  by 


V  -  V  COS (Y„  -  O  ) 


V 

^itch 


V  cos  (Y  ) 

”  ’“pitch 


(3.23) 


The  acceleration  components  are  then  a  f\inction  of  the  angular 
rate  of  change  of  the  velocity  vector 


a  -  V  Y 

^pitch  ^itch 


Y 

®y*ir  ®yair 


(3.24) 


The  angular  acceleration  commands  are  then  distributed  to  the 
missile's  cartesian  coordinate  accelerations  using  the 
following  geometric  relationships 
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x'  = -a  sin  (a  cos  (a  ) 

‘^iteh  "Utah  pitch'  yaw' 

y‘  = -a  sin  (a  _  . )  sin  (a  ) 

Vtch  •^Itoh  Pitch^  yaw' 


z  =  a  cos  (a  A 

•Utcb  "Wtcb  pitch' 


(3.25) 


X  =  -a  sin  (a  ) 

®y.»,  V*' 

y  =  a  cos (o  ) 

“y«v 


The  missile  acceleration  in  each  plane  is  then 


X  =  X 


+  X 


y  =  y  +  y 

*"  ”V>ltch  “yw 


(3.26) 


z  =  z 


and  the  total  missile  acceleration  is 


««  ■  ♦  y'a  * 


(3.27) 


2 .  Commeuid  Guidance 

The  CLOS  autopilot  also  takes  the  guidance  commands  and 
translates  them  into  missile  accelerations.  Similar  to  the 
proportional  navigation  autopilot,  this  autopilot  translates 
the  angular  accelerations  into  cartesian  coordinate 
accelerations . 

The  commanded  angular  accelerations  of  equation  (3.15) 
are  translated  to  cartesian  accelerations  using  the  following 
relationships 
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^n,  = s±n(a  Jcosio  ) 

”^ltcb  <=^lteh  pitch'  yaw' 

yL  = sin  (a  sin  (a  ) 

"Wteh  ‘=^lteh  ^  pitch'  yaw' 

=  a  .  cos  (a 

"Wtefc  <=^iee6  pitch' 

X  =  -a  sin  (a  ) 

andy^  yaw' 


(3.28) 


V  =  a  cos  ( o  ) 


The  overall  missile  accelerations  are  also  given  by  equations 
(3.26)  and  (3.27)  . 


F.  MISSILE  AMD  TARGET  KINEMATICS 

The  missile  and  target  kinematics  can  be  developed  using 
the  state  space  representation 


x=fx  X  y  y  z  zl’ 

jn  [jn  m  ^  m  ^  m  m  mj 

X.  .  [x.  y.  y.  z,  f  J- 


(3.29) 


The  system  is  then  represented  by 


X  =  A  X  +  B  u 

/n  in  tn  m  nt 


Xj.  =  A^  Xj.  +  B^ 


(3.30) 
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where 


0  1  0  0  0  0 
0  0  0  0  0  0 
0  0  0  1  0  0 
<^"  0  0  0  0  0  0 
0  0  0  0  0  1 
0  0  0  0  0  0 


(3.31) 


A  signal  flow  graph  for  the  missile  and  target  kinematics  can 
be  seen  in  Figure  11. 


G.  KALMAN  FILTER  DEVELOPMENT 

The  introduction  of  noise  into  the  simulation  creates  a 
more  realistic  scenario.  The  problem  is  to  determine  the 
target's  flight  path  by  filtering  the  noise.  This  simulation 
uses  an  extended  Kalman  filter  to  estimate  the  target's 
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flight . 

The  noisy  observed  target  flight  is  the  input  to  the 
filter.  The  cartesian  and  spherical  coordinates  of  the  target 
are  then  used  in  the  Kalman  iteration  to  estimate  the  target's 
position.  The  filter  is  developed  to  use  preprocessed  linear 
pseudomeasurements.  These  measurements  are  given  by 

x{kT) 

y{kT) 


where 

a  =  LOS  pitch  angle 
3  =  LOS  yaw  angle 

The  measurement  equation  then  becomes 

1  0  0  0  0  0 

y{kT)  =  0  0  1  0  0  0  x{kT)  *  (3.33) 

0  0  0  0  1  oj 

where  =  N(0,R),  and  R  =  H  (kT)  R'h’' (kT)  .  H(kT)  and  R*  are 
given  by 
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H(kT)  = 


5x(r,a, 3) 

hx(r,ct,  3) 

5x(r,a, 3) 

5r 

5a 

53 

6y(r,a,3) 

5y(r,a, 3) 

5y(r,a, 3) 

5r 

5a 

53 

5z(r,a, 3) 

bz{r,ct,  3) 

5z(r,a, 3) 

6r 

5a 

53 

R- 


a  2  0  0 

r 

0  a  2  0 

a 

0  0 


(3.34) 


The  discrete  time  system  model  then  becomes 


x(  (Jt+1)  T)  =  Fx(Jtr)  + 


«  N(0,Q) 


Q  = 


E  = 


ZOO 
0  Z  0 
0  0  Z 

gT^  qT' 
“1“ 

qT^ 


2 

qT 


(3.35) 


The  initial  condition  for  the  filter  is 


x(0)  «  W(x;,Po) 


(3.36) 


The  Kalman  algorithm  is  then  given  by  the  following  set  of 
equations 


K. 


x". , 

Jt-l/Jt  k/k 

=  FP^.^F^  *  Q 

k*l/k  k/k 

.  K^R^K/ 


(3.37) 
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IV.  SIMULATION  RESULTS 


A.  OVERVIEW 

The  proportional  navigation  and  CLOS  simulations  are 
tested  using  three  target  flight  scenarios.  In  the  first 
scenario  the  target  has  constant  velocity  and  level  flight  in 
two  dimensions.  In  the  second,  the  target  has  a  constant 
velocity  and  level  flight  in  three  dimensions.  Finally,  in 
the  third,  noise  is  added  to  the  three-dimensional  scenario. 

The  Simulink  models  and  associated  MATLAB  code  for  the 
proportional  navigation  and  CLOS  simulations  are  contained  in 
the  Appendix. 

B.  SIMULATION  ASSUMPTIONS 

The  following  assumptions  are  held  throughout  the 
simulation : 

(1)  Acceleration  due  to  gravity  does  not  effect  the 
missile  or  the  target. 

(2)  The  missile  is  lying  in  the  xy  plane  at  launch. 

(3)  Missile  acceleration  is  limited  to  30g. 

(4)  The  proportional  navigation  constant  is  N=6 . 
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C.  SIMULATION  SCENARIOS 


1 .  Constemt  Velocity  In  Two  Dimensions 

The  first  scenario  is  a  two-dimensional  engagement.  The 
target  is  flying  at  a  constant  altitude  with  no  acceleration. 
The  target  parameters  are  as  follows 


30000  ft 

-3000  ft/s 

^'t 

0  ft/s  2 

yt 

0  ft 

yt 

0  ft/s 

y't 

0  ft/s  2 

Zt 

1000  ft 

0  ft/s 

z’t 

0  ft/s^ 

2 .  Constemt  Velocity  In  Three  Dimensions 

The  next  scenario  is  a  three-dimensional  engagement.  The 
target  is  flying  at  a  constant  altitude  with  no  acceleration. 
The  target  parameters  are  as  follows 

-  60000  ft 
}t^  -  -2121  ft/s 

-  0  ft/s^ 

-  10000  ft 

-  -2121  ft/s 
.  0  ft/s  2 

Zj  -  1000  ft 

Zj.  -  0  ft/s 

-  0  ft/s^ 
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3.  Three-Dimensional  Simulation  With  Radar  Noise 

The  final  simulation  uses  the  same  target  parameters  as 
the  three-dimensional  constant  velocity  simulation.  White 
noise  is  added  to  the  target  flight.  This  simulates  received 
noise  in  the  target's  radar  return.  The  noise  has  the 
following  characteristics 


-  15  ft 


yaw 


pitch 


D.  RESULTS  AND  SIMULATION  COMPARISON 

Figure  12  indicates  the  missile  leads  the  target.  This 
is  attributed  to  the  slow  missile  autopilot  time  constant 
(T;y,=l  sec)  and  the  target's  speed  advantage  of  mach  3  to  mach 
2  over  the  missile.  This  problem  is  exaggerated  in  figures  12 
and  13  since  the  z  scale  is  twenty  times  the  x  scale.  It  was 
determined  by  considering  the  z  acceleration  profile  in  figure 
16,  the  z  velocity  profile  in  figure  19,  and  the  z  position 
profile  in  figure  12,  that  this  effect  was  caused  by  the 
autopilot . 

Figure  14  shows  the  rate  of  change  of  o  is  positive  for 
approximately  1  second;  thereafter  it  is  negative  but,  for  2 
seconds  the  missile  has  a  positive  commanded  acceleration. 
Figures  15  and  16  show  the  missile's  acceleration  variations. 
Figures  17,  18  and  19  show  the  missile's  velocity  variations. 
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Target  Trajectories  in  the  xz  Plane. 


Dimensional  Engagement  Plot. 


28 


0  1  2  3  4  5  6 

Time  (seconds) 


Figure  14.  Proportional  Navigation  Scenario  1.  LOS  Angle 
Opitch . 


Figure  15.  Proportional  Navigation  Scenario  1.  Commanded 
Acceleration  in  the  x  Direction. 
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Figure  16.  Proportional  Navigation  Scenario  1.  Commanded 
Acceleration  in  the  z  Direction. 


Figure  17.  Proportional  Navigation  Scenario  1.  Total  Missile 
Velocity. 


Figure  18.  Proportional  Navigation  Scenario  1.  Missile 
Velocity  in  the  x  Direction. 


Figure  19.  Proportional  Navigation  Scenario  1.  Missile 
Velocity  in  the  z  Direction. 
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Plots  for  the  other  scenarios  are  given  in  the  Appendix. 
The  following  table  summarizes  the  missile's  closest  point  of 
approach  (CPA) ,  and  the  time  of  the  CPA  for  each  simulation. 


Scenario 

Simulation 

CPA 

Time  of 

CPA 

1 

Prop  Nav 

4.13  ft 

5.89 

s 

CLOS 

1.39  ft 

7 .18 

s 

2 

Prop  Nav 

14.94  ft 

14  .72 

s 

CLOS 

1.24  ft 

19 . 51 

s 

3 

Prop  Nav 

27.15  ft 

14 .5 

s 

CLOS 

267.79  ft 

22.34 

s 

Table  1 .  Missile  Miss  Distance  Summary 


Overall,  the  proportional  navigation  missile  achieves  a 
quicker  target  intercept  time.  The  miss  distances  for  each 
missile  are  very  close,  except  when  noise  is  added.  The  CLOS 
missile  degrades  significantly  in  the  presence  of  noise. 

The  proportional  navigation  missile  is  a  superior 
missile.  The  CLOS  missile  is  unable  to  give  satisfactory 
results  when  sensor  noise  is  added  to  the  simulation.  For 
very  short  range  intercept  scenarios,  where  sensor  noise  is 
negligible,  the  missile  will  perform  well.  The  proportional 
navigation  missile  will  perform  well  for  any  engagement 
scenario.  This  fact  makes  proportional  navigation  preferable 
for  missile  guidance. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 
A.  CONCLUSIONS 

The  simulation  provides  insight  in  chosing  the  proper 
type  of  missile  guidance.  The  two  types  of  guidance  explored 
both  give  acceptable  miss  distances  without  sensor  noise 
However,  when  sensor  noise  is  present  the  proportional 
navigation  missile  outperformed  the  CLOS  missile. 

The  presence  of  an  on  board  seeker  gives  the  proportional 
navigation  missile  an  advantage  when  dealing  with  sensor 
noise.  Since  the  sensor  is  on  the  missile  as  it  closes  the 
target,  the  sensor  noise  will  have  less  of  an  effect  on  the 
detection  of  the  target.  The  CLOS  missile  is  guided  from  a 
stationary  radar  at  the  launch  site.  The  error  incurred  from 
sensor  noise  does  not  decrease  as  the  missile  approaches  the 
target.  To  overcome  this  problem  the  CLOS  missile  will 

require  a  very  sophisticated  tracking  radar  that  has  very 
little  sensor  noise. 

The  addition  of  noise  to  the  engagement  provides  a  more 
realistic  scenario  for  the  missile  control  problem. 
Developing  a  noise  filter  and  adjusting  the  missile 
characteristics  to  adapt  to  the  noise  created  a  unique  and 
educational  challenge.  The  increased  realism  reinforced  the 
fact  that  actual  missile  control  developement  is  a  compromise 
of  design  requirements. 

B .  RECOMMENDATIONS 


The  simulation  can  be  taken  to  several  different  levels. 
The  target  flight  can  be  modified  for  different  engagement 
scenarios.  A  manuevering  target  would  provide  another  level 
of  realism  to  the  engagement. 

An  adjoint  model  could  be  built  for  each  simulation. 


33 


This  would  aid  in  the  miss  distance  analysis  for  the  two 
missiles. 

Finally,  different  noise  filters  can  be  developed  and 
tested.  The  miss  distance  will  be  decreased  if  better  noise 
filtering  is  achieved  during  the  simulation. 
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APPENDIX 


A.  COMMAND  GUIDED  MISSILE  PLOTS  FOR  SCENARIO  1 


Missile/Target  Engagement  in  the  xz  Plane 


Figure  20.  Command  Guidance  Scenario  1.  Missile  and  Target 
Trajectories  in  the  xz  Plane. 
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Plot. 


Figure  22.  Command  Guidance  Scenario  1.  Cross  Range  Error. 
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Figure  23.  Command  Guidance  Scenario  1.  Commanded 
Acceleration  in  the  x  Direction. 


Acceleration  in  the  z  Direction. 


3 


Figure  25.  Command  Guidance  Scenario  1.  Total  Missile 
Velocity. 


Figure  26.  Command  Guidance  Scenario  1.  Missile  Velocity  in 
the  X  Direction. 
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fl/sec 


the  z  Direction. 
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Distance  (ft) 


B.  PROPORTIONAL  NAVIGATION  MISSILE  PLOTS  FOR  SCENARIO  2 


Proportional  Navigation  Missile/Target  Engagement  Scenerio  2. 


Distance  (ft  x  1  e4)  “2  0  Distance  (ft  x  1  e4) 


Figure  28.  Proportional  Navigation  Scenario  2.  Missile  and 
Target  Trajectories  in  Three  Dimensions. 
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Figure  29.  Proportional  Navigation  Scenario  2.  Missile 
Target  Trajectories  in  the  xy  Plane. 
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Figure  31.  Proportional  Navigation  Scenario  2.  Missile 


Target  Trajectories  in  the  yz  Plane. 


Figure  32.  Proportional  Navigation  Scenario  2.  LOS  Angle 
Opitch . 
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Figure  33.  Proportional  Navigation  Scenario  2.  LOS  Angle 
Oyaw. 


Figure  34.  Proportional  Navigation  Scenario  2.  Commanded 
Acceleration  in  the  x  Direction. 


3 


3  5  10  15 

Time  (seconds) 

Figure  35.  Proportional  Navigation  Scenario  2.  Commanded 
Acceleration  in  the  y  Direction. 


Commanded  Acceleration  z  Direction 


Figure  36.  Proportional  Navigation  Scenario  2.  Commanded 
Acceleration  in  the  z  Direction. 
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fl/sec 


5 


0 


Missile  Velocity  in  the  y  Direction 


Figure  39.  Proportional  Navigation  Scenario  2.  Missile 
Velocity  in  the  y  Direction. 


Missile  Velocity  in  the  z  Direction 


Figure  40.  Proportional  Navigation  Scenario  2.  Missile 
Velocity  in  the  z  Direction. 
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Distance  (ft) 


C.  COMMAND  GUIDED  MISSILE  PLOTS  FOR  SCENARIO  2 


Command  To  Line  Of  Sight  Missile/Target  Engagement  Scenerio  2 


Distance  (fix  1e4)  0  Distance  (fix  1e4) 


Figure  41.  Command  Guidance  Scenario  2.  Missile  Target 
Engagement . 
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2  3  4  5  6 

Distance  (ft  x  1  e4)  ^  q4 


Figure  43.  Command  Guidance  Scenario  2.  Missile  Target 
Trajectories  in  the  xz  Plane. 
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Figure  44.  Command  Guidance  Scenario  2.  Missile  and  Target 
Trajectories  in  the  yz  Plane. 


Figure  46.  Command  Guidance  Scenario  2.  Cross  Range  Error  in 
the  Pitch  Plane. 


Figure  47.  Command  Guidance  Scenario  2.  Cross  Range  Error  in 
the  Yaw  Plane. 
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Commanded  Acceleration  in  the  x  Direction 


Figure  48.  Command  Guidance  Scenario  2.  Commanded 
Acceleration  in  the  x  Direction. 
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Figure  50.  Command  Guidance  Scenario  2.  Commanded 
Acceleration  in  the  z  Direction. 


Figure  51.  Command  Guidance  Scenario  2.  Total  Missile 
Velocity. 
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53 


Figure  54.  Coinmand  To  Line  Sight  Scenario  2.  Missile 
Velocity  in  the  z  Direction. 
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Distance  (ft) 


PROPORTIONAL  NAVIGATION  MISSILE  PLOTS  FOR  SCENARIO  3 

Missile  Target  Engagement 

1500 

1000 

500 

0 
2 

xIO^  -2 

Distance  (ft  x1e4)  ^  0  Distance  (ft  x1e4) 

Figure  55.  Proportional  Navigation  Scenario  3.  Missile  and 
Actual  Target  Trajectory. 
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Figure  56.  Proportional  Navigation  Scenario  3.  Missile  and 
Filtered  Target  Trajectory. 


Distance  (ft  x  1e4) 


Figure  57.  Proportional  Navigation  Scenario  3.  Missile  and 
Target  Trajectories  in  the  xz  Plane. 
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0  1  2  3  4  5  6  7 

Distance  (ft  x  1  e4)  ^  o'* 


Figure  58.  Proportional  Navigation  Scenario  3.  Missile  and 
Target  Trajectories  in  the  xy  Plane. 


Figure  59.  Proportional  Navigation  Scenario  3.  Missile  and 
Target  Trajectories  in  the  yz  Plane. 
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Figure  61,  Proportional  Navigation  Scenario  3.  Missile 
Accelerations . 
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E.  COMMAND  GUIDED  MISSILE  PLOTS  FOR  SCENARIO  3 


Trajectory. 
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0  1  2  3  4  5  6  7 

_  Distance  (ft  x  1  e4)  „  ^ 


Figure  63.  Command  Guidance  Scenario  3.  Missile  and  Target 
Trajectories  in  the  xy  Plane. 


Distance  (ftx1e4) 


Figure  64.  Command  Guidance  Scenario  3,  Missile  and  Target 
Trajectories  in  the  xz  Plane. 
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Trajectories  in  the  yz  Plane. 


I _ _ _ Time  (sec) _ 

Figure  66.  Conimand  Guidance  Scenario  3.  Cross  Range  Error. 
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F.  PROPORTIONAL  NAVIGATION  SIMDLINK  MODEL 
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Figure  68.  Target  Flight  with  Noise. 
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Figure  69.  Target  Flight  without  noise. 


Figure  70.  Proportional  Navigation  Pitch  and  Yaw  Seeker. 
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Figure  71.  Proportional  Navigation  Pitch  and  Yaw  Autopilot. 


sin (u) 


Figure  72.  Proportional  Navigation  Missile  Flight  in  the  x 
Direction. 
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Figure  73.  Proportional  Navigation  Missile  Flight  in  the  y 
Direction. 


Figure  74.  Proportional  Navigation  Missile  Flight  in  the  z 
Direction. 
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Figure  75.  Proportional  Navigation  CPA  and  Stop  Simulation 
Calculation. 


Figure  76.  Proportional  Navigation  CPA  Calculation. 
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G.  COMMAND  GUIDANCE  SIMULINK  MODEL 


Figure  77.  Command  Guidance  Missile  Model. 
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Figure  78.  Command  Guidance  Target  Flight  with  Noise. 
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Figure  80.  Command  Guidance  Missile  Flight. 
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target  sqrt  (xm^2+ym''2+zin""2)  Rt 

flight 


Figure  81.  Command  Guidance  Target  Range. 


Flight 


Figure  82.  Command  Guidance  Sigma  Yaw  Calculation. 
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flight 


Figure  83.  Command  Guidance  Sigma  Pitch  Calculation. 


Figure  84.  Command  Guidance  Sigma  Calculation. 
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Figure  85.  Command  Guidance  Cross  Range  Error. 


Figure  86.  Command  Guidance  Cross  Range  Error  Yaw. 
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crep 


sxgma 

pitch  target 
sigma 

pitch  missile 


MATLAB 

Function 


sign(u[l]’-u[2]  ) 


CRE%  ='Jrt(u[l)^2-u(2)^2) 

- ►du/dt - 

>|crepd| 

Derivative 

CREdot  pitch 

Figure  87.  Command  Guidance  Cross  Range  Error  Pitch. 
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Figure  89.  Cartesian  to  Spherical  Block  for  both 
Simulations. 
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H.  MISCELLANEOUS  MATLAB  CODE 

^  program  generates  the  noise  used  in  the  target  flight. 

randn{ 'seed' ,26579) ; 

ti=  [0  ;  .001:30]  ; 

for  i=  1:30001 

%Range  noise 

U(i)  =  randn*15; 

%Pitch  angle  noise 
V(i)  =  randn*pi/l80; 

%Yaw  angle  Noise 

W(i)  =  randn*pi/l80; 
end 

%This  program  sets  the  initial  conditions  for  the  Kalman 

%Filter .  It  is  run  at  the  beginning  of  each  simulation. 

clear  P 

clear  xhat 

global  P 

global  xhat 

%initial  covariance  matrix 
P=le6*eye  (6 )  ; 

%initial  estimated  target  position 
xhat=  [10000  -500  1000  -500  0  500]'; 

%This  function  runs  a  Kalman  filter  algorithm 

%for  the  given  A,  B,  C  matices  for  constant  velocity  flight 

function [xhat , P] =klmn(u, P,xhat)  ; 

%initiali2ation 
A=[0  1  0  0  0  0; 

000000; 

000100; 

000000; 

0  0  0  0  0  1; 
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0  0  0  0  0  0]; 

B=[0;0;0;0;0;0]; 

C=[1  0  0  0  0  0; 

0  0  1  0  0  0  ; 

0  0  0  0  1  0]; 

%Time  step  and  q  parameter  for  the  Kalman  Filter 
q=l; 

dt=  .001; 

%Continuous  to  discrete  conversion 
[phi, del] =c2d(A,B,dt) ; 

%Specify  position  and  angle  vectors 
pos=  [u(l)  ;u(2)  ;u(3)  ]  ; 
ang=  [u(4)  ;u(5)  ;u{6)  ]  ; 

%Calculate  Sigma  Matrix 

sigma= [ ( (q*2) * (dt) *3) /3  ((q‘2)*(dt)*2)/2; 

( (q‘2)* (dt) *2)/2  (q*2)*(dt)]; 

%Calculate  Q  Matrix 

0= [sigma  zeros (2, 2)  zeros (2, 2)  ; 
zeros(2,2)  sigma  zeros(2,2); 
zeros{2,2)  zeros(2,2)  sigma]; 

%Calculate  Rstar  Matrix 

Rstar=  [225  0  0;0  {I*pi/180)  *2  0;0  0  (l*pi/l80)  *2] 
%Kalman  iteration 
xhatkpl  =  phi*xhat; 

Pkpl  =  phi*P*phi'  +  Q; 

H  =  calch(ang) ; 

R  =  H*Rstar*H'; 

K  =  Pkpl*C' *inv(C*Pkpl*C' +R)  ; 

Pklkl  =  (eye(6) -K*C) *Pkpl* (eye(6) -K*C) '  +  K*R*K' ; 

xklkl  =  xhatkpl  +  K* (pos ( : , 1) - (C*xhatkpl) ) ; 
global  xhat ; 
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xhat  =  xklkl; 
global  P; 

P  =  Pklkl; 

%This  function  calculates  the  H  matrix  for  computation 
%in  a  Kalman  filter 
function  [H]=calch(x) 
a= (tan(x(2) ) ) “2; 
b= ( tan(x (3) ) ) *2 ; 
c= (sec (x (2) ) ) *2 ; 
d= (sec (x (3) ) ) *2 ; 
e=tan (x(2)); 
f =tan(x (3) ) ; 

H  =  [1/ ( (a*b+a+b+l) * . 5)  -x (1) * (b*e*c+c*e) / ( (a*b+a+b+l) *1 . 5) 

-X (1) * (a*d*f +d*f ) / ( (a*b+a+b+l) “1.5) ; 

e  /  (  (a*b  +  a  +  b  +  l)  *  .5) 

-X  ( 1)  *  (a*c*  (b+1)  )  /  (  (a*b+a+b+l)  *1 . 5)  +x  (l)  *c/  (  (a*b+a+b+l)  “  . 5) 
-X (1) *e*f *d* (a+1) / ( (a*b+a+b+l) *1.5)  ; 

f/((l+b)*.5)  0  -x(l)*b*d/((l+b) “1.5)+x(l)*d/((l+b) *.5)]  ; 
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